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ABSTRACT 



Aims. The occurrence of new FU Orionis-like objects is fundamental to understand the outburst mechanism in young stars and their 
role in star formation and disk evolution. Our work is aimed at investigating the properties of the recent outburst of V1647 Ori. 
Methods. Using optical and mid infrared long slit spectroscopy we monitored VI 647 Ori in outburst between February 2004 and 
O |. January 2006. 

Results. The optical spectrum is characterized by Ha and H/3 in P-Cygni profile and by many weak Fel and Fell emission lines. 
Short timescale variability was measured in the continuum and line emission. On January 2006 we detected for the first time forbidden 
emission lines ([OI], [S II] and [Fell]). These lines are likely produced by an Herbig-Haro object driven by V1647 Ori. The mid 
infrared the spectrum of V1647 Ori is flat and featureless at all epochs. The SED changed drastically: the source was much redder in 
the early outburst than in the final phase. The magnitude rise and the SED of VI 647 Ori resembles that of a FUor while the duration 
and recurrence of the outburst resemble that of a EXor. The optical spectrum is clearly distinct from either the absorption line spectrum 
of a FUor or the T Tauri-like spectrum of an EXor. 

Conclusions. Our data are consistent with a disk instability event which led to an increase of the mass accretion rate. The data also 
suggest the presence of a circumstellar envelope around the star+disk system. The peculiar N band spectrum might be explained by 
dust sublimation in the outer layers of the disk. The presence of the envelope and the outburst statistics suggest that these instability 
events occur only in a specific stage of a Class I source (e.g. in the transition phase to an optically visible star surrounded by a 
£^ , protoplanetary disk). We discuss the outburst mechanisms in term of the thermal instability model. 

\^») , Key words. Protoplanetary disks - Accretion Disks - Instability - Herbig Haro object 
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' 1. Introduction complex, showed an increase of its optical/IR brightness of up 

| to 6 magnitudes. The outburst has been observed from the X-ray 

Li , One of the clearest pieces of evidence for disk accretion during regime (e g Qr osso et al. l2005T l to infrared wavelength (see e.g. 

. early stages of stellar evolu tion ar e events like FU Ono nis and M uzerolle et al. 120051 Andrews et al. l2004l . In February 2004, 4 

»H I EX Lupi outbursts (Herbig [19771 Hartmann & Kenyon 1996). months after the onset of me outburst , the brightness rise stopped 

«? , These outbursts are thought to be the consequence of a sudden and the magn i tude remained (relatively) constant. In November 

and steep increase of the mass accretion rate onto the central star, 2 005, a communication by Kospal et al. d2005T l claimed the be- 

which changes from those commonly found around T Taun stars ginning of a fast fading phase in the optical light of vl647 Qri 

(-10 M Q yr ) into values of 10 -10 M G yr .Statistical The system is further characterized by a red energy distribution 

studies suggest that young low-mass stars experience several FU and by many emission lines in its optical and near-IR spectrum. 

Ononis like outbursts during the early phase of stellar evolution. Apart from the Bracket t series (seen purely in emission), all the 

The emergence of a new pre-main-sequence outburst object is Hydrogen lines exhibit a P-Cygni profi le, which indicates mass 

then a unique opportunity to address the physical process that outflow in a wind Vacca et al (2004) find that their nea r-IR 

occurs in the disk's interior. For this reason several astronomers emission fine spectrum is consistent with a dense and ionized 

have focused their attention on the recent outburst of V1647 On. wind mode i where the optically thick H lines are produced. The 

V 1 647 Ori is a young eruptive star known to be the illuminat- sam e model is able to explain why t he outb urst has not been seen 

ing source of McNeil's Nebula, a reflecting nebula discovered by at ra dio wavelength (Andrews et al. l2004b . V1647 Ori is known 

the a mateu r astronomer J.W. McNeil in January 2004 (McNeil to have experienced a previous outburst as is clear from the ap- 

et al. 2004). In the months following the discovery, the star, lo- pearance of the reflection nebula in 1996-1997 in the atlas of 

cated in the LI 630 cloud within the Orion B molecular cloud |Mallas & Kreimer (1970)| and as recently confirmed by Aspin et 

al. ( 120061 1. Furthermore, its optical and near-IR spectrum does 

Send offprint requests to: D. Fedele not resemble any other previous spectra of FUors or EXors ob- 

* Based on observations collected at the European Southern jects (Reiburth & Aspin 2004, Vacca et al.|2004). The 2-3 years 

Observatory, Paranal, Chile. Proposal ID: 272.C-5045, 074.C-0679, duration of the outburst, its recunence on a timescale of decades 
276.C-5022 
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and the "peculiar" spectrum of VI 647 Ori, are important clues 
for the comprehension of outburst events in pre-main-sequence 
stars. 

From February 2004 to January 2006 we have followed the 
evolution of the outburst of VI 647 Ori at optical (4700-7300 
A) and mid-IR (8-13 fim) wavelengths. Here, we present the 
results of our photometric and spectroscopic monitoring of the 
eruptive source. The analysis of the reflection nebula are pre- 
sented in a separate paper (Fedele et al. 2007 hereafter paper 
I). Observations and data reduction are described in section [2] 
In section [3] we analyze the observations. A discussion of the 
results is presented in section|4] We draw the conclusions in[5] 

2. Observations and data reduction 

Observations were performed using FORS2 at ESO's Very Large 
Telescope in Paranal, Chile and TIMMI2 at the 3.6 m tele- 
scope at La Silla. FORS2 ( Appe nzeller et al. 1998] > is an opti- 
cal facility (3000-10000 A) which allows imaging in different 
bands and grism spectroscopy. TIMMI2 (Kaufl et al. 2003]) is a 
mid-infrared (8-14 fim) multi-mode instrument including low- 
and medium-resolution spectrograph. We have also included in 
our analysis some publicly available VLTI/MIDI observations 
of VI 647 Ori taken from the ESO archivefl MIDI is the mid- 
infrared beam-combiner facilities of the ESO VLT interferome- 
ter (Leinert et al. |20031 >. 

2. 1 . Optical Spectroscopy 

12 long slit spectra of VI 647 Ori were obtained with FORS2 
and with the grism 1400V (4560-5860A, A/AA ~ 2100) between 
2004 February 18 and 2005 December 27. Seven further long 
slit spectra were obtained between 2004 December 08 and 2006 
January 29, with the grism 1200R (5750-7310A, A/AA ~ 2100). 
A detailed log of the observations is reported in Table Q] 

A standard optical long slit spectra extraction procedure was 
applied to reduce the raw data - bias subtractions, flat-fielding, 
wavelength calibration, cosmic rays and sky background re- 
moval and weighted average along the spatial axis. Observations 
of spectro-photometric standard stars during each night allowed 
us to compute the sensitivity function of the spectrograph. To 
flux-calibrate the spectra of VI 647 Ori, the 1 -dimensional ex- 
tracted spectra were first multiplied by the sensitivity function 
- to compute the exact slope of the spectrum - and then scaled 
to the flux level measured from the acquisition images. This is 
straightforward for the red spectra since the acquisition images 
were taken with the same filter. For the blue spectra we have as- 
sumed (V-R c ) ~ 1 .8+0.2 (as measured by Mc Gehee et al. 120041 
in February- April 2004 and from Kospal et al. 2005 in October 
2005) from which we computed the absolute V flux level. 

Differential aperture photometry of VI 647 Ori has been 
computed over an aperture radius of 2'.' 5 2 (10 pixels) from the 
acquisition images (Rc filter). In paper I we computed the Rc 
magnitude of two references stars in the FORS2 frames which 
were used to calibrate the instrumental magnitude of the acqui- 
sition images discussed here. The two stars were found to be not 
variable and to have the following Rc magnitude: RA(J2000) 
= 05:46:09.71; DEC(J2000) = -00:03:31.1; R c = 20.08 + 0.09 
andRA(J2000) = 05:46:05.88; DEC(J2000) = -00:02:39.7; R c = 
16.39 ± 0.02. The latter corresponds to the comparison star "E" 
in Semkov [2003 and |2006l who measured R c = 16.39 ± 0.02. 
The results are listed in Table Q] 



2.2. Mid-IR spectroscopy 
2.2.1. TIMMI2 

Low-resolution (A/AA a 200) /V'-band (7.7-13.0 /an) spectra of 
V1647 Ori (IRAS 05436-0007) were obtained on 2004 March 
08 (JD 2453072.526) and 2006 January 10 (JD 2453746.663) 
using the TIMMI2 instrument on the ESO 3.6 m telescope at La 
Silla. Sky subtraction was achieved by chopping in the North di- 
rection with an amplitude of 10" followed by a nodding pattern 
with the opposite direction and amplitude. The total integration 
time was 25 minutes per spectrum. The slit, with a width of 172, 
was centered on a compact source detected in the A^'-band acqui- 
sition image. A spectrum of the reference star HD 37160 (K0III) 
was obtained before or after each IRAS 05436-0007 observa- 
tion. Data were reduced using the usual steps of residual back- 
ground subtraction, spectral extraction, and wavelength calibra- 
tion. Correction for the telluric ozone absorption bands, as well 
as absolute flux calibration were achieved by ratio-ing the IRAS 
05436-0007 spectrum to that of HD 37160 flux-calibrated us - 
ing the spectral templates by Cohen et al. (Co hen et al. 1 999). 
The resulting spectra are shown in Fig. [6] 



2.2.2. VLTI/MIDI 

We also analyzed interferometric observations of VI 647 Ori ob- 
tained with VLTI/MIDI on 3 nights between 2004 December 30 
and 2005 March 01. The interferometric data were previously 
analyzed by Abraham et al. (2006). Hence we will discuss here 
only the MIDI low-resolution (R=30) spectrum acquired on each 
night by the instrument after the interferometric observation. 

Using the MIA+EWS-1.30 software package, a fixed mask 
was applied to the MIDI chopped spectrum. The background has 
been estimated from the off-source (sky) frames and then sub- 
tracted from the on-source frames. Finally the one-dimensional 
spectrum has been extracted. In the same way, spectra of MIDI 
calibrator stars were extracted. 

Aperture photometry was computed from the acquisition im- 
ages with a narrow filter centered at 8.7 /jm adopting a fixed 
aperture of 12 pixel (l'.'O). The throughput of the two MIDI 
channels are, for an unknown reason, different, and the aperture 
photometry differs for the two telescopes. We used the results 
from the channel B which tend to be the more stable of the two. 
Conversion factors from counts to Jy were computed from the 
MIDI calibrators, whose 8.7 /mi flux were evaluated from the- 
oretical spectral energy distributions. Energy distributions have 
been derived from the Cohen list (Coh en et al. 19 99 ) of infrared 
standard stars for the calibrator HD 37160 (F^j^, = 11.5 Jy), 
and by matching the spectral type of the calibrator to stars in 
the Cohen list and scaling the spectrum with the ratio of their 
IRAS 12 micron fluxes for HD 107446 (F 8 7/Jm = 37.9 Jy) and 
HD 50778 (Fsj^m = 28.9 Jy). The three calibrators were also 
used to compute the sensitivity function of MIDI (wavelength 
dependence of instrument's response). The spectrum of V1647 
Ori has been flux calibrated by multiplying it by the sensitivity 
function and by scaling it in order to match the 8.7 fim flux. The 
result is shown in Figure [6] 



1 http://www.eso.org/archive 



- http: //www. strw. leidenuniv. nl~koehler/MIA +EWS- Manual 
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JD -2 450 000 


FWHM 


Slit width 


Slit PA 


Spectral range 


Exposure time 


S/N 


Rc 


(UT) 




(") 


(") 


o 


(A) 


(sec) 




(mag) 


2004-02-18 


3053.081 


0.85 


1.0 


45.0 


4560-5860 


700 


20 


17.39 + 0.10 


2004-02-23 


3058.090 


1.25 


1.0 


45.0 


4560-5860 


700 


19 


17.31 ±0.10 


2004-03-13 


3077.018 


0.9 


1.0 


45.0 


4560-5860 


700 


14 


17.52 ±0.10 


2004-03-18 


3082.011 


1.15 


1.0 


45.0 


4560-5860 


700 


13 


17.42 ±0.10 


2004-03-27 


3091.020 


0.9 


1.0 


45.0 


4560-5860 


700 


14 


17.31 ±0.10 


2004-12-08 


3347.315 


1.2 


0.7 


90.0 


4560-5860 


500 


9 


17.23 ± 0.05 


2004-12-21 


3360.251 


0.65 


0.7 


90.0 


4560-5860 


500 


14 


16.87 ± 0.05 


2005-01-05 


3375.180 


1.35 


0.7 


90.0 


4560-5860 


500 


10 


17.07 ± 0.05 


2005-02-18 


3419.126 


1.0 


0.7 


90.0 


4560-5860 


500 


4 


17.77 ± 0.05 


2005-02-29 


3430.114 


0.65 


0.7 


90.0 


4560-5860 


500 


9 


17.12 ±0.05 


2005-03-15 


3444.045 


0.72 


0.7 


90.0 


4560-5860 


500 


7 


17.50 ± 0.05 
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a v i cnn 
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ZU. /4 ± U. 1 1 


2004-12-08 


3347.322 


1.15 


0.7 


90.0 


5750-7310 


500 


20 


17.23 ± 0.05 


2004-12-21 


3360.258 


0.68 


0.7 


90.0 


5750-7310 


500 


42 


16.87 ± 0.05 


2005-01-05 


3375.187 


1.08 


0.7 


90.0 


5750-7310 


500 


26 


17.07 ± 0.05 


2005-02-18 


3419.133 


1.12 


0.7 


90.0 


5750-7310 


500 


15 


17.77 ± 0.05 


2005-02-29 


3430.121 


0.62 


0.7 


90.0 


5750-7310 


500 


23 


17.12 ±0.05 


2005-03-15 


3444.052 


0.69 


0.7 


90.0 


5750-7310 


500 


22 


17.50 ± 0.05 


2006-01-29 


3764.122 


0.73 


0.7 


90.0 


5750-7310 


4 x 1800 


4 


22.05 ±0.11 



3. Results 

3.1. Optical lightcurve 

Figure[TJshows the light curve of V 1 647 Ori in the Rc band based 
on the data of Table Q] and on previous measurements by other 
authors. Due to the influence of McNeil's nebula on the com- 
putation of the stellar flux and subtraction of local background, 
measurements with different instruments may result in a differ- 
ent magnitude estimation. In particular, given the better spatial 
resolution of our data, we may better disentangle the contribu- 
tion from VI 647 Ori from that of the nebula. For this reason, an 
offset of -0.3 mag was applied to Rc measurements by other au- 
thors. Such offset is not needed for the data of Acosta-Pulido et 
al. 2007. Three further Rc measurements from paper I are plot- 
ted in Figure [TJ 

Since measurements of Rc of the early outburst are not avail- 
able in literature, we estimated Rc from Ic measurements of 
Briceno et al. 2004 The (Rc - Ic) color, as measured by var- 
ious authors (McGehee et al. 2004, Ojha et al. 120061 1. seems rea- 
sonably stable during the plateau and the fading phase, showing 
(Rc - Ic) ~2-0 with a scatter of 0.2 mag. With this we calcu- 
late the expected Rc magnitudes for the period October 2003 - 
February 2004 from the Ic measurements, assuming the same 
(Rc - Ic) color for the rising phase. From figure [TJ the opti- 
cal light curve of V1647 Ori can be divided in three parts: i) 
from October 2003 to February 2004 - the rising period; ii) from 
February 2004 to August 2005 - the plateau phase and iii) from 
August 2005 to January 2006 - the fading period. 

The rising part is very steep: from October 2003 to January 
2004 the optical magnitude increases by more than 3 magnitudes 
in Rc- From the pre-outburst magnitude level, Rc ~ 23.5, com- 
puted by McGehee et al (see inset in figure [1} it results that the 
total rise in brightness of VI 647 Ori is > 6 mag in Rc- From 
the light curve in figure [TJ we find a rate of increase of Rc of 
~ 1.5 mag/month. Assuming that this rate remained constant 
during all the rising phase, we estimate that the outburst began 
around the middle of August 2003, slightly earlier than October- 
November 2003 as found by Briceno et al. (2004). 

During the plateau phase the optical brightness shows a slow 
decline with time (ARc = 0.02mag/month), on top of which Rc 



displays a non-periodic, flickered, oscillation on short timescale. 
The light curve in figure QJ combines data from different works 
which have used different instruments and aperture sizes. This 
may result in a systematic scatter of the data from work to work. 
However, an intrinsic variation of the optical brightness on short 
time scale is clearly present. From our data we measure a vari- 
ation of ~ 0.5 mag between 2004 Feb. 17 (R c = 16.90 ± 0.05, 
paper I) and Feb. 18 (R c = 17.39 + 0.10, Tableffl). Thus, V1647 
Ori at its maximum light shows an optical brightness variation on 
a time scale of 24 hours. For five nights we have two consecutive 
acquisition images (separated by a few minutes) from which we 
searched for even smaller time scale variations of Rc, however, 
no significant change in optical brightness (ARc > 0.1 mag) are 
detectable from these measurements. The short timescale vari- 
ability is of the same order of that found by other authors in 
optical and near infrared (Ojha et al. 120061 Walter et al. 12004b 
on a timescale of a week. The total duration of the plateau phase 
is less than 2 years. 

Our photometry confirms the rapid brightness decrease an- 
nounced by Kospal et al. (2005). Four months after the claimed 
onset of the fading phase (dated to August 2005, when Orion 
reappeared on the sky, however from the light curve it is clear 
that it started before), the brightness of VI 647 Ori was still di- 
minishing. From August 2005 to January 2006 Rc dropped of 
4 mag. On 2006 January 29, the last Rc measurements, we es- 
timate Rc = 22.05 + 0.11 which is still more than 1 magnitude 
above the pre-outburst level of McGehee et al. ((2004). From the 
light curve we estimate a fading rate of ~ 0.8 mag/month dur- 
ing this phase. Assuming a constant fading rate, Rc reached the 
pre-outburst level at the beginning of April 2006. 

It is worth to note that the light curve is not symmetric, the 
rising and fading phase have a different slope. From the two dif- 
ferent rates we may infer that the physical process which led 
to the beginning of the outburst is not simply the reverse of 
the process ending this outburst. The rate of brightness varia- 
tion during the three phases are in agreement with that found by 
|Acosta-Pulido et al. (2007)1 
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Fig. 1. Optical, Rc, light curve of V1647 Ori. Data are from: filled circles this work; squares Semkov, E.H. (2006 2004); triangles 
Kospal et al. ( 120051); star s |Acosta-Pulido et al. (2007)[ dia monds Briceno et al. ( 120041 1 ; filled squares McGehee et al. 020041 ); open 
circles Ojha et al. (2006}; star |Acosta-Pulido et al. (2007)[ An offset of -0.3 mag was applied to the imported data (apart from those 
of Acosta-Pulido et al.2007) in order to reach a similar magnitude level with our data. In the inset the pre-outburst magnitude level 
reported by McGehee et al. (2004) is shown. 



3.2. Optical spectra 

The positive slope of the optical spectrum of V1647 Ori (fig.|2]i 
reveals a red energy distribution of the source. Clearly visible 
are the Ha and H/3 lines both characterized by a P-Cygni profile. 
The He I ,15875, Nal Dl & D2 doublet in absorption and Fel 
and Fe II lines in emission (P-Cygni) are detected in the spectra 
taking during the plateau phase. Due to bad pixel columns in 
correspondence of [OI] A5517 and [O 1] /16300 lines a residual 
of the data reduction is present in the spectra taken during the 
plateau phase. For this reason we cannot find evidence of any 
source features at these wavelengths. 

During the plateau phase the optical spectrum shows only 
minor changes on equivalent width and line flux (see discus- 
sion below). The overall shape of the spectrum remains "con- 
stant" during this period. The spectrum taken during the fad- 
ing phase (blue - 2005 December 27, red - 2006 January 29, 
Fig.|2]i is completely different from all the others: superimposed 
on a very faint continuum with a slightly positive slope there 
are very strong emission lines in correspondence of Ha, [O I] 
X16300,6363, [S II] ^6717,6731 and [Fell] at 7172 A . The 
sky line [OI] A5511 is saturated and it is not possible to find 
evidence of source emission at this wavelength. 

From the flux-calibrated spectra, for each line detected, we 
computed the equivalent width (EW) and the line flux (F/) mul- 
tiplying EW by the continuum level at the line center (Tables [3] 
&|4|. A major source of uncertainties in this computation is the 



determination of the continuum level, which in turn depends on 
the accuracy of the sensitivity function and of the aperture pho- 
tometry from the acquisition images. The final accuracy is of the 
order of 10% on EW and of 15% on F,. 

The P-Cygni profile of Ha and H/5 is in both cases asym- 
metric with the emission components lacking the high velocity 
part (see Figures [3j |4|. This profile results from deviation of 
spherically-symmetric wind and is observed in FU Ori objects 
and T Tauri stars and can be explained by the presence of an 
opaque disk which occults part of the redshifted emission (see 
e.g. Hartmann 1998). The profiles of the two lines differs sig- 
nificantly (see figures [3] & 0): the H/J has a strong and wide ab- 
sorption and a "weak" narrow emission while the Ha has a huge 
emission and a weak absorption. In both cases, the blue-shifted 
absorption shows at least two components: one at — 450 km s _1 , 
and the other at — 150 km s . While the low velocity compo- 
nent remains almost constant over all the plateau phase, the high 
velocity one is variable. In particular, the latter shows in both 
lines a progressively decrease in extension from February 2004 
to March 2005 until the whole absorption disappear in the fading 
phase spectrum. Furthermore, on three nights (2005 January 05, 
February 18 and March 15) the bluest absorption component of 
the Ha is "replaced" by an emission. Also the emission compo- 
nent varies from night to night, displaying a change in equivalent 
width and line flux. 

P Cygni signatures are also displayed by Fe lines. However, 
due to the low S/N of the spectra, the absorption component is 
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Fig. 2. Examples of optical spectrum of VI 647 Ori obtained during the plateau phase (top two panels), and during the fading phase 
(bottom two panels). Plateau spectrum: clearly visible are the Ha and Wf3 with P-Cygni profile, Fel (5328, 6191, 6495 A) and 
Fell (5169, 6432, 6516 A) lines and absorption from Nal Dl & D2 (5889, 5895 A) and Hel line (5875 A). A residual of the data 
reduction is present in correspondence of [OI] (6300, 5577 A) due to bad pixel columns, and between 6870 - 6910 A due to a 
non-perfect removal of telluric absorption. Fading phase spectrum: no lines are detected in the blue part ([OI] 5577 A is a residual 
of the data reduction), while the red part is characterized by a strong emission from He, [O I] 6300, 6363 A, [S II] 67 1 7, 673 1 A and 
[Fell] 7172 Alines. 



clearly detected only for the Fe II A5 169.08 transition. The max- 
imum absorption is at — 370 km s _1 and has a terminal velocity 
of ~ - 600 km s _1 . The other Fe emission lines show variation 
in equivalent width and line flux with time. 

Thus, similar to what we have found for the continuum emis- 
sion (R c magnitude), the emission lines vary on time scales of 
months and weeks. If continuum and line emission are produced 
in the same region, then the flux of a generic emission line (Fi) is 
linearly proportional to continuum flux (F R - R c bandpass flux). 
In figures [5] we plotted the line flux of the strongest emission 
lines detected (Ha and H/3) versus F«. Arrows indicate upper 
limits. Assuming a power law dependence of the line emission 
from F« (Fi oc FT) we searched for a correlation between line 
and continuum emission. The continuous line in figures [5] is the 
best fit to the data. We find ju a — 0.8 and yfip — 0.7. The Ha 
best fit is more robust than the H/3 one (the Ha emission is much 
stronger and it is also detected during the fading phase). All the 
other emission lines detected are weak and strongly affected by 
the low S/N of the spectra. The power law model fit in this case 
is more uncertain and clear evidence of dependence of the line 
emission from the continuum emission cannot be found. The 
value of y found for the H Balmer lines is close to unity. This 
means that the variations in line flux is correlated to the vari- 



ation of the continuum emission. This evidence might suggest 
that the two emissions arise from the same region of the system, 
or at least, that the emission mechanisms are physically linked. 

Contrary to what was found by |Acosta-Pulido et al. (2007)] 
we do not find any clear trend of the Balmer emission lines with 
time during the plateau phase. Both the EW and the line flux 
show a random variation with time in this period. 

For the first time we detected optical forbidden lines in the 
spectrum of VI 647 Ori taken during the fading phase (see figure 
12 bottom right). The spectrum clearly shows strong emission 
from [OI] /Li6300,6363, [S II] Ai6717,6731, [Fell] A7172 an 
Ha. These lines provide evidence for hot (a few thousand K) gas 
close to V1647 Ori. These are used as tracers of Herbig-Haro 
objects where a collimated jet from the central star collides with 
the ambient medium. The emission is produced by the cooling of 
the shocked gas. It is interesting to note that Eisloffel & Mundt 
(fl997T > have already identified IRAS05436-0007 as the driving 
source of HH 23. Similar and perhaps newly formed ejecta could 
be responsible for the forbidden emission lines seen here. None 
of these forbidden lines were previously detected in the plateau 
spectrum, most likely because of the overwhelming continuum. 



D. Fedele et al.: Temporal evolution of the outburst of VI 647 Ori 











2004-12-21 




?no5 m os 




2005-02-18 


2.5e-17 






2.5e-17 






2.5e-17 




2.5e-17 




1.9e-17 






1.9e-17 






1.9e-17 




1.9e-17 






1.3e-17 






1.3e-17 






1.3e-17 


I 


1.3e-17 






6.3e-18 






6.3e-18 






6.3e-18 


VJ 


6.3e-18 






Oe+00 




0e+00 






0e+00 




0e+00 






i — H ' 


-1000 


1000 


-1000 1000 


-1000 1000 


-1000 C 


1000 


's 

i 




















2005-02-29 






2005-03-15 




2006-01-29 







1.5e-17 






1.5e-17 


l.lc-17 






l.le-17 


7.5e-18 






7.5e-18 


3.8e-18 






3.8e-18 


0.0e+00 




0.0e+00 




2.5e-19 


L 




1.9e-19 








1.3e-19 








6.3e-20 








0.0e+00 




-1000 


(1 


1000 



v [km s 



-1, 



Fig. 3. Time evolution of the Ha emission in the period December 2004 - January 2006. In all the "plateau spectra" the line has an 
asymmetric P-Cygni profile. In the "fading phase spectrum", the line is pure in emission with no trace of absorption. 



Table 2. Mid-infrared photometry of VI 647 Ori from TIMMI2 
and MIDI data. For comparison in column 3 and 4 we report 
the corresponding Rc magnitude and the Rc-Nl color, where 
Nl is the magnitude at 8jum converted to the ESO mid infrared 
photometric system (I Van der Bliek et al. 19961 1. 



DATE 




Rc 


Rc-Nl 


Ref. 




[Jy] 


[mag] 


[mag] 




2004-03-07 


6.4 


16.8 


14.4 


Muzerolle et aU2005b 


2004-03-08 


6.5 


16.8 


14.4 


this work (TIMMI2) 


2004-03-11 


4.3 


16.8 


14.0 


Andrews et al. d2004l> 


2004-12-31 


2.3 


16.7 


13.2 


this work (MIDI) 


2005-01-03 


2.5 


16.7 


13.3 


this work (MIDI) 


2005-02-21 


2.3 


17.1 


13.5 


this work (MIDI) 


2005-03-01 


2.5 


17.2 


13.7 


this work (MIDI) 


2006-01-10 


0.5 


21.0 


15.8 


this work (TIMMI2) 



3.3. Mid-Infrared spectra 

The rise in brightness of VI 647 Ori during the outburst has 
been witnessed by others authors also at longer wavelength (see 
e.g. A ndrews et al. 1.2004., Muzerolle et al. 120051 Abraham et al. 
2006) . Our TIMMI2 spectra (see fig. [6]) confirm the increased 
mid-infrared flux: from the pre-outburst level of 0.53 Jy at 12 
/mi (IRAS, Point Source Catalog) up to 7.6 Jy on 2004 March 
08. The 8 - 14/mi spectrum is essentially featureless and flat 
all along the spectral range. This result is not consistent with 
the mid-infrared spectrum in Andrews et al. (2004) taken with 
UKIRT/Michelle only three days after our TIMMI2 spectrum. 
Their measurement reveals a strong red energy distribution with 



the flux going from 4.0 Jy at 8 pm up to 12.0 Jy at 12p.m. Our 
estimate of the mid-infrared flux is however consistent with the 
SPITZER/IRAC observations of Muzerolle et al. J2005l > taken 
on 2004 March 07. They measure indeed a flux of 6.43 Jy at 8 
pm, very close to our estimate. A residual of the data reduction 
in correspondence of the strong atmospheric absorption bands 
centered at 9.6 (O3) and 12.55 /mi (CO2) do not allow us to bet- 
ter analyze these two regions. 

Ten months later, on December 2004, the mid-infrared flux 
of VI 647 Ori revealed by the MIDI observations dropped by a 
few Jy. The spectrum is again flat and featureless. The MIDI 
spectra analyzed here, were all taken during the optical plateau 
phase (Fig. |6j. They all reveal a flat and featureless spectrum. 
Within the accuracy of these spectra (10%), the flux level remain 
constant between December 2004 and March 2005. Thus, also in 
the mid-IR the system experienced a plateau phase. 

The rapid optical brightness fading is also experienced by the 
system in the mid-infrared: on 2006 January 11, the flux level 
of our TIMMI2 spectrum at 12 pm is 0.9 Jy, still considerably 
higher than the pre-outburst level. Also in this case the spectrum 
is flat and featureless. 

It is worth to note that the optical and mid-IR light curve 
of VI 647 Ori are different. In the optical the brightness in- 
creases during the rising phase remaining below the value of the 
plateau phase and finally decreases. On the other hand, the mid- 
IR brightness is higher during the "mid-IR rising phase" than 
during the "mid-IR plateau phase" (figure|7]l. 
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Fig. 4. Time evolution of the H/3 emission in the period February 2004 - December 2005. As for the Ha, during the plateau phase 
the line has a P-Cygni profile. No H/? emission is detected in the low S/N "fading phase spectrum". 
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Fig. 5. Line flux of Ha and HJ3 versus continuum flux. The best power law fit is over-plotted. 

4. Discussion 



Pre-main-sequence stars are known to be intrinsically variable 
objects. The variability mechanisms might be different: solar- 
like coronal activity, spots on the stellar surface, stellar pulsa- 
tion, partial obscuration and subsequent clearing of the line of 
sight. These processes are however unable to generate the ~ 44 
L luminosity increment produced by V 1 647 Ori and to produce 



the > 6 optical magnitude burst lasting more than 2 years. To 
release such an amount of energy the existence of a secondary 
luminosity source is necessary. Similar brightening events from 
FU Orionis stars are explained by a sudden increase of the mass 
accretion rate from a circumstellar disk onto the central star. 
The increased accretion rate produces an accretion luminosity 
(L aC c M) which may overwhelm the stellar brightness. Such a 
process can explain both the dramatic brightening (from X-ray 
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Table 3. Lines detected in the blue spectra of V1647 Ori between February 2004 and December 2005. EW is negative for emission 
lines and positive for absorption lines. In case of emission lines also the line flux is reported EW are expressed in A and line flux in 
10~ 18 W irT 2 . The accuracy on equivalent width and line flux is respectively of the order of 10% and 15%. 



Went. 


A[A] 
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2004-02-23 


2004-03-13 


2004-03-18 






EW 




EW 


f line 


EW 


^Une 


EW 


F/,„ c 


FV?(em) 


4861 


-2.91 


1.6 


-1.94 


1.34 


-5.38 


3.12 


-2.22 


1.26 


H/?(ab) 




5.45 




5.91 




11.48 




5.68 




Fe 1 1 (em) 


5169 


-0.58 


0.51 


-0.17 


0.18 


-1.68 


1.52 


-0.37 


0.32 
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0.74 




1.04 
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Fig. 6. Time evolution of the mid-Infrared spectrum of VI 647 
Ori taken with TIMMI2 (black and blue lines) and MIDI (red 
line). The mid-infrared spectrum measured by MIDI remain con- 
stant, within the uncertainties (10%), between December 2004 
and March 2005. Here, only the MIDI spectrum taken on 2005 
February 21 is shown. 

to the infrared) as well as the strong Ha emission observed in 
the recent outburst of V1647 Ori. Kastner et al. (2006) recently 
indeed confirmed that the X-ray evolution of VI 647 Ori in out- 
burst reflects the near-infrared evolution and is consistent with 
the hypothesis of an increased mass accretion rate. 

As a consequence of the enhanced accretion rate a strong 
wind develop from the disk's surface. The blue-shifted absorp- 
tion component of Ho- and Hf3 in the spectrum of VI 647 Ori 



Fig. 7. 8 fim light curve of V1647 Ori. Filled triangles - 
TIMMI2 data; filled squares - MIDI data; open circle - 
SPITZER/IRAC data from Muzerolle et al. d2"0"05T l; open trian- 
gle - UKIRT/MICHELLE data from Andrews et al. ([20041 . The 
inset shows the pre-outburst IRAS measurements at 12 fj.m. 

are signatures of this wind. The disappearance of the absorption 
component in the Ha during the fading phase is a confirmation 
that the strong wind ceased and that the system has been going 
back to a phase of slow accretion. In this system the origin of the 
H Balmer emission lines is controversial since both the wind and 
the mass infall might contribute to the line formation. The mag- 
netospheric accretion predicts a further contribution to the lines 
emissivity produced in the optically thin free-falling accretion 
columns. 
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Table 4. Continuation of table[3]for red spectra between December 2004 and January 2006. 
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A[A] 


2004-12-08 


2004-12-21 


2005-01-05 


2005-02-18 






EW 




EW 


p 


EW 


p 


EW 


p 


He I 


5875 


1.14 




0.71 




0.32 




<1.6 




Nal (Dl) 


5889 


1.44 




0.97 




2.26 




2.82 




Nal (D2) 


5895 


1.01 




0.70 




1.17 




2.88 




Fel(169) 


6191 


-0.31 


0.54 


-0.35 


-0.74 


-0.34 


0.63 


>0.08 


<0.8 


[01] 


6300 


















[OI] 


6363 














- 


- 


Fell 


6432 


-0.74 


1.58 


-0.44 


1.17 


-0.45 


1.03 


-0.48 


0.58 


re 1 


o4bO 


-1.08 


2.46 


-0.55 


1.54 


-0.73 


1.76 


-0.83 


1.08 


Cp TT 


DJ ID 


-0.55 


1.27 


-0.43 


1.24 


-0.53 


1.30 


-0.69 


0.93 


Ho'(em) 


OJDZ 


-26.61 


66.34 


-28.74 


88.66 


-30.27 


80.44 


-38.81 


56.60 


rlQ^aD ; 




3.30 




2.55 




0.74 




<0.53 




TQ TT1 

La llj 


fn 1 n 
o / 1 / 














- 


- 


tti 


O / J 1 














- 


- 


rp TTl 

[be llj 


1 1 1 L 














- 


- 


Went. 


A[A] 


2005-02-29 


2005-03-15 


2006-01-29 










EW 


F/,„ e 


EW 




EW 








He I 


5875 


<0.8 
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0.87 




1.51 




<7.3 








Fel(169) 


6191 


-0.24 


0.44 


-0.30 
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-40.98 


0.42 
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0.08 






T- TT 

re 11 
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-0.42 


0.94 
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0.85 
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<0.03 






Fel 
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-0.52 


1.23 


-0.92 


1.59 


>-2.7 


<0.03 






Fell 


6516 


-0.35 


0.85 


-0.52 


0.92 


>-2.7 


<0.03 






Hff(em) 


6562 


-17.25 


44.65 


-22.07 


39.94 


-108.30 


1.66 






Ha(ab) 




1.78 




1.10 




<3.9 








[SII] 
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-4.52 


0.09 






[SII] 
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-7.24 


0.15 






[Fe II] 
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-9.48 


0.32 







The Fell lines at 6432.68 and 6516.81 A detected in the 
spectra of VI 647 Ori both originate from multiplet 40. They are 
commonly also seen in emission in the spectra of strongly ac- 
creting young stars such as V380 Ori (Rossi 1999; Shevchenko 
[19991 , Z CMa (Hessman et al. [19911 Garcia et al. [19991 van den 
Anck er et al. 120041 , PV Cep and MWC 1080 (Hernandez et al. 
2004) and are amongst the strongest emission lines in such en- 
vironments. These Fe II lines appear to be absent from the spec- 
trum of young stars showing lower accretion rates. The presence 
of Fel 6495.81 emission and the absence of strong Fell lines 
from higher multiplets suggests that the iron-line forming region 
in VI 647 Ori may be somewhat cooler than the typical electron 
temperatures of 10,000-20,000 K inferred in for example Z CMa 
(van den Ancker et al. |2004| l. 

4.1. The circumstellar envelope and the puzzling mid 
infrared spectrum 

The accretion disk alone is not able to produce the long wave- 
length (A > 10 fim) emission observed, unless it flares strongly 
over a large range of distance scale (see e.g. Hartmann 1998). 
The sub-millimeter continuum flux during the outburst remains 
at its pre-outburst level and there are no signatures of flux 
changes in these wavelength regime (Andrews et al. (2004)). 
These findings are consistent with the presence of a dusty cir- 
cumstellar envelope, probably a remnant of the infalling en- 
velope. Intriguingly, Kenyon & Hartmann ( 119911 ) suggested 
the presence of the infalling envelope to explain the enhanced 
mid-infrared flux from FU Orionis objects in eruptive phase. 
Muzerolle et al. (120051 hereafter M04) attempt to reproduce the 



spectral energy distribution (SED) of VI 647 Ori by means of a 
standard viscous accretion disk and of an optically thin envelope. 
Their model predicts a 10/im emission feature that is produced 
by silicate dust grains. However, our multi-epoch mid-infrared 
spectroscopy reveals a flat and featureless spectrum during the 
whole outburst duration (see Fig. [6]). This is highly unusual. In 
FU Orionis objects the silicate feature is seen sometimes in emis- 
sion (VI 057 Cyg, FU Ori, BBW 76, V1515 Cyg, Green et al. 
2006) and sometimes in absorption (V346 Nor, Z CMA Green 
et al. 120061 Acke & van den Ancker 2004). These differences 
are probably caused by differences in the optical thickness of the 
system (disk + envelope) at 10 fim. 

The model proposed by Abraham et al. (2006) adopts a 
simple viscous accretion disk model (without envelope). Their 
model predicts a flat and featureless mid infrared spectrum. 
However, the emission at longer wavelength (> 10/im) requires 
an highly flared disk. In such a model the outer part of the disk 
is directly illuminated during the outburst producing a flux en- 
hancement also in the (sub-)millimeter. This is not observed 
(Andrews et al. 2004 ). The stability of the (sub-)millimeter emis- 
sion favors the presence of dusty circumstellar envelope. 

The excitation plot for the CO fundamental ro-vibrational 
lines (Rettig et al. 120051 ) and the detection of Av = 2 CO band 
heads at ~ 2.3 fxm (Vacca et al. |2p04) indicate the presence of 
hot (~ 2500 K) and dense gas. Since the dust sublimes at ~ 1500 
K, such emission likely arises in regions free of dust. Our sug- 
gestion is that even in the mid infrared the bulk of the emis- 
sion is produced by the gas in a dust-free region of the disk. 
Nevertheless we cannot exclude a contribution from the dust (e.g 
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Table 5. Typical value of outburst from pre-main-sequence stars. The outburst recurrence of FU Orionis objects has been estimated 
has the time needed to replenish the disk mass after an outburst with a constant infall rate. 





FUors 


EXors 


VI 647 Ori 


Outburst duration [yr] 


> 10 


~1 


2.6 


Outburst recurrence [yr] 


>200 


5-10 


37 


Mass accreted during an outburst[M Q ] 


> 10~ 3 


io- 6 - 10~ 5 


2.5- 10 s 


Magnitude variation [optical mag] 


4-6 


2-5 


~6 


Accretion Luminosity [L Q ] 


few 10 2 


>25 


44 


Outburst accretion rate [M Q yr -1 ] 


10~ 4 


10~ 6 - 10~ 5 


10~ 5 


Envelope infall rate [M Q yr -1 ] 


5-10~ 6 


10~ 7 - 10~ 6 


7-10~ 7 


Wind velocity [Km s -1 ] 


>300 


200-400 


300 - 400 


Mass loss rate [M Q yr -1 ] 


10~ 6 - 10~ 5 


10~ 8 - 10~ 6 


4-10" 8 


Spectral features 


absorption spectrum 


emission line spectrum, 


emission line spectrum, 




F/G-type supergiant like 


T Tauri like, Ho- inverse P Cyg 


strong Ha emission (P Cygni) 




deep CO absorption 


CO abs./em., Bry emission 


CO abs./em., Bry emission 
Forbidden lines in fading phase 


Note 






X-rays variability 



arising at larger radii in the disk). The emission at longer wave- 
length is dominated by the dust in the envelope. 

Figure [7] shows the temporal evolution of the mid infrared 
flux during the outburst. The flux decay at these wavelengths 
is faster than the decay in the optical (compare with figure Q]). 
The 8-14 jjm flux drastically decreases from March to December 
2004 while the optical continuum remains constant over the 
same period. This produces radical changes in the SED of V1647 
Ori during the outburst. As can be seen from the Rc - Nl color 
(Tabled the system is redder during the early outburst. A likely 
explanation is that in the earlier phases of the outburst the disk 
was hotter and a larger region of the disk contributed to the emis- 
sion seen at 10 /mi. However, a detailed modeling of the obser- 
vations is necessary to explain which parameters (such as tem- 
perature, opacity, surface area) are responsible for the observed 
changes. 

5. Conclusion 

Outbursts in pre-main-sequence stars have been historically clas- 
sified in two main groups upon their similarity to the prototypes 
FU Orionis and EX Lupi (Herbig 1 19771 1 depending on outburst 
duration, maximum magnitude variation and spectral features at 
maximum light. Table[5]lists the main characteristics of the two 
groups and that of VI 647 Ori. 

Since the onset of the outburst of VI 647 Ori it has been de- 
bated whether this system is either a FUor or an EXor object. 
VI 647 Ori resembles some aspects of an EXor (outburst dura- 
tion, recurrence of the outburst), and some aspects of a FUor 
(magnitude rise, SED). However the recurrence timescale of 
the outburst has intermediate value between the two classes. Its 
emission line spectrum is clearly distinct from either the absorp- 
tion line spectrum of a FUor or the T Tauri-like spectrum of an 
EXor (where the H lines show inverse P-Cygni profile). 

VI 647 Orionis is not the only outburst PMSs suspected of 
having an intermediate nature between the two main classes. OO 
Ser experienced recently an outburst which lasted ~ 5-10 years 
(Kospal et al. 2006), too fast for a FUor and too slow for an 
EXor. The SED of OO Ser is typical of a FU Orionis object, and 
has roughly the same shape in quiescent and outburst phase. 

A common denominator in all young eruptive stars detected 
so far seems to be the presence of circumstellar material as well 
as that of a reflection nebula. These structures are likely rem- 
nants of the inf ailing envelope. The infalling envelope is a po- 
tential reservoir of mass for the disk which experiences recursive 



outbursts. Assuming that Lj, / during the outburst is dominated 
by the accretion luminosity, M04 estimate a mass accretion rate 
of ~ 10~ 5 M yr~'. Considering the 2-3 years duration of the out- 
burst and the 37 years recurrence timescale, a constant envelope 
infall rate of ~ 7-10~ 7 M yr~' is necessary to replenish the disk 
after each outburst. The disk accretion rate during the quiescent 
phase is estimated to be ~ 6-10~ 7 M yr~' (see e.g. M04). 

Submillimeter maps reveal that FU Orionis stars have ac- 
cretion disks that are larger and more massive than those of T 
Tauri stars (Sandell & Weintraub 2001) and are comparable in 
mass to those seen around Class I sources (i.e. young stellar ob- 
jetcs with flat or rising infrared SED). The circumstellar material 
ar ound V I 647 Ori accounts for 0.04 + 0.01 M (Tsukagoshi et 
al. |2005) which is slightly larger the disk mass of a T Tauri star 
(~ 0.01 M ). All these findings suggest that outbursts occur in 
Class I sources, where the star is still embedded in the infalling 
envelope. The outburst duration and mass accretion rate during 
outburst seem to correlate with the infall rate (see Table [5}: ob- 
jects with higher infall rate have longer outburst and reach higher 
accretion rate while objects with smaller infall rate experience 
short-lived outbursts. The occurrence of short outbursts might 
suggest that the envelope is becoming more and more thin, i.e. 
that the system is in a transition phase of an embedded Class I 
source to an optically visible star surrounded by a protoplanetary 
disk (Class II). 

The Orion Nebula Cluster (ONC) is subject to extensive ob- 
servational campaigns. So far, roughly 1600 stars have been con- 
firmed to be members of the ONC and 55% of these (at least) 
posses a circumstellar disk (Hillenbrand et al. 1998). Assuming 
that Class I sources account for 20% of stars with disks (lower 
limit, equal to the Class I fraction found in the original paper of 
Lada & Wilking 1984), we expect to have more than 200 Class I 
sources in the ONC only. If all of these 200 sources experience 
FU Orionis-like activity, we would expect to see 10-40 (de- 
pending on outburst duration and recurrence of FUors, VI 647 
Ori and EXors) of these stars in outburst at any time in the ONC. 
The total number of young eruptive stars discovered so far in 
the whole Orion star forming region is only seven: three FUors 
(FU Orionis itself, V883 Ori and Reipurth 50), three EXors (NY 
Ori, VI 1 18 Ori, VI 143 Ori) and V1647 Ori. Not all of these are 
in outburst at the same time. It appears that there is a deficit of 
observed outbursts in Orion. We conclude that not all Class I 
sources undergo FU Orionis-like events for their entire lifetime. 

A possible solution to this problem is that we are over- 
estimating the number of "true" Class I sources in Orion. The 
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spectral energy distribution of an isolated T Tauri source seen 
with a disk close to edge-on may mimic a typical Class I SED. 
An alternative solution is that outbursts occur only in a specific 
stage of the early evolution, namely, in the transition phase of 
an embedded Class I source to an optically visible T Tauri or 
Herbig AeBe star. 

We also caution that this result is based only on a small por- 
tion of Class I sources in Orion, namely, those not embedded 
in high density regions. Recent Spitzer observations show that 
there is a large number of such sources which he in regions of 
extremely high extinction (Megeath et al. 2006 ). An outburst in 
one of these objects could easily have been missed in the existing 
surveys of Orion. More regular infrared surveys of star forming 
regions are badly needed to investigate the number of FUor-like 
outbursts in these embedded sources. 

5. 1 . Outburst mechanism 

Instability mechanisms of different flavors have been proposed 
to explain the FU Orionis phenomenon. The gravitational forces 
of a companion star may perturb the disk enhancing accretion 
(see e.g. Bonnell & Bastien 1992). Gravitational instability has 
been also proposed but, for this instability to occur the disk has 
to be cold and massive. The most accepted mechanism to trigger 
such outbursts is via thermal instability in the inner disk (see 
e.g. Lin & Papa loizou [19851 Clar ke et al. [19901 Kawazoe & 
Mineshige 119931 Bell et al. 119951 ). A key point of this model 
is that of an high accretion rate in the outer disk, of the order 
of few x 10~ 6 M yr _1 . According to this model, outbursts will 
occur as long as mass is deposited in the outer disk at such high 
rate. This implies that the outbursts will become shorter in time 
and smaller in amplitude as the infall into the disk ceases. An al- 
ternative explanation for eruptive events has been suggested by 
Gammie (1996). He suggested that "dead zones" of decreased 
accretion may develop in the case of magnetic viscosity. Material 
from the outer part of the disk may accumulate at the edge of 
such dead zones until high-accretion-rate episodes occur. 

One way to explain the different outburst properties of 
FUors, EXors and VI 647 Ori is by introducing a different out- 
burst mechanism for each of these classes of objects. However, 
the data collected so far in all young eruptive stars raise the pos- 
sibility of the existence of only a unique class of outburst objects. 
In this scenario, the different types of outburst are produced by 
a continuum variation of one or more of the parameters involved 
in the instability, rather than a variation of physical mechanisms. 
For instance, the presence of a remnant of the infalling envelope 
and the estimated infall rate (see Table [5]) might favor the ther- 
mal instability as the unique model to explain the three groups. 
If this is the case, the different outburst duration between a FU 
Orionis, VI 647 Ori and an EX Lupi might be explained with 
the difference in the infall rate. The thermal instability model 
predicts indeed shorter and smaller amplitude outbursts as in- 
fall ceases below 10~ 6 M yr~' to approach typical T Tauri disk 
accretion rate of ~ 10~ 7 M yr _1 , 

The same model is able to explain also the diversity in spec- 
tral features of the three classes of objects. The mass accretion 
rate reached during the outburst is (also) dependent on the en- 
velope infall rate. When very high accretion rates and disk tem- 
peratures are reached, as in the case of FUors, the disk's interior 
may become hotter than the disk surface and the emitted spec- 
trum will show absorption lines. If the accretion rate is not high 
enough to invert the temperature gradient in the disk's interior, 
the final spectrum will be dominated by emission features which 
are produced in the disk's wind and/or in the magnetospheric 



accretion columns. In this case we will see EXors or V1647 Ori- 
like spectra. 
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